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ABSTRACT 


Presently available viral dynamic research models are of the non-linear mixed effect models Most of the non-linear models 
developed by using the differential equations. Finding the solution of parameters are difficult if it is non-linear model. So, 
In this paper an attempt has been made to derived expected time of new viral budding consider kolmogrov in the life cycle 
of viral Dynamics. Predictive distribution of viral replication using modified Exponential power distribution through the 


Bayesian methodology has been derived. 
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1. INTRODUCTION 


At the end of 2011, globally, 34.0 million people were living with HIV. An estimated 0.8% of adults aged 15-49 year are 
living with HIV. The world Health organization (WHO) states clinical failure indicating severe immune deficiency among 
adults and adolescents. HIV produces millions of its copies in the bloodstream and consequently weakens the immune 
system by defeating WBCs. This type of WBC is called CD}T cells or T-helper cells. HIV infected persons must go 
through antiretroviral treatment to fight infections. There is a way to treat HIV and no Vaccine is available yet, but 
HAART is an effective technique used to slow down the progression of the disease. The immunological failure as CD} T 
count falls to the baseline (below 100 cell/mm’) virology failure as plasma viral load measurements after 3 month with 
adherence support.CD}T Count is best predictor for disease status and immediate risk of death and also used to identify those 
who have advanced HIV disease.There is no medicine to cure HIV infection is a challenge of physicians. So the researcher 
has developed statistical and mathematical modelling for the viral prediction. There are challenges in statistically modelling 
immune responses to longitudinal HIV viral load exposure as function of covariates. Bayesian Markov chain Monte Carlo 
mixed effect models to incorporate priors and examine different distributional assumptions. The improvement of health in 
HIV/AIDS patients on Highly Active Antiretroviral the reply is characterised by an increase in CD{T counts and decrease in 
viral load to undetectable levels. In this research, seven stages of viral life cycle considered as continues time, discrete stage 
markov chain and modelling of viral replication in the Bayesian approach. In modelling HIV/AIDS progression in patients, 
researchers mostly deal with either viral load or CDT cell counts only. The orthogonal covariates viral load along with 


CDj{T baseline and using Bayesian methodology, viral prediction is determined. 
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2. REVIEW OF LITERATURE 


Onoja Matthew Akpa (2010)studied the some of the models proposed by various authors for describing the epidemiology 
as well as the epidemiological consequences of the HIV/AIDS epidemic and how some of them could be modified to suit 
the situations in other countries. They also discussed the limitations and the place of such models in the fight against the 
HIV epidemic.R. Kannan and M. lIyappan (2017)have studied a stochastic model for estimating the expected time to 
seroconversion of HIV infected using preventive strategy. The effectiveness of the adoption preventive strategy in the 
elongation of the expected time to seroconversion by introducing the stochastic model under the assumption that the 
threshold level of antigenic diversity as a random variable which follows exponentiated modified Weibull distribution. G. 
Meenakshi, S. Lakshmi Priya (2018)have used the Bayesian Methodology is used to estimate the HIV replication for the 
future period per infected cell. In the HIV life cycle, there are seven different stages, such as binding, fusion, reverse 
transcription, integration, replication, assembly and budding. At the stage of binding, the various viral components bind 
with the cell. Such as a-Surface glycoprotein (gp120),b - trans membrane glycoprotein (gp41) CCRS5, c - reverse 
transcriptase (P66 /P51) are considered to estimate the HIV replication. New risk function derived to use the components a, 
b and c. Finally, posterior parameter of HIV replication in a cell is estimated. The estimated parameter is compared with 
some standard parameter using various loss functions.G. Meenakshi and S. Saranya (2018) have derived a new modelled 
estimate the HIVreplication periodically using Markov processes in the condition of decay of CD{Tcells. The proposed 
model is illustrated through numerically. Ali Raza, Muhammad Rafiq(2019) has explained a reliable numerical analysis of 
a stochastic HIV/AIDS model in a two-sex population considering counselling and antiretroviral therapy (ART). The 
authors are comparing the solutions of the stochastic and deterministic HIV/AIDS epidemic model. Xiaodan Sun, Hiroshi 
Nishiura(2019)has discussed Estimating HIV incidence is crucial for monitoring the epidemiology of this infection, 
planning screening and intervention campaigns, and evaluating the effectiveness of control measures. However, owing to 
the long and variable period from HIV infection to the development of AIDS and the introduction of highly active 
antiretroviral therapy, accurate incidence estimation remains a major challenge. Claris Shoko and Delson Chikobvu (2019) 
have studied CD}T count and viral load as an essential component in monitoring HIV treatment outcome. However, CD{T 
cell count monitoring sometimes fails to predict virological failure resulting in unnecessary switch of treatment lines which 
causes drug resistance and limitations of treatment options. Yan Wang, Tingting Zhao (2019) stochastic perturbation is 
introduced into HIV model with virus-to-cell infection, cell-to-cell transmission, CTLimmune response and _ three 
distributed delays. The stochastic integro-delay differential equation is transformed into degenerate stochastic differential 
equations. Examine that the influence of random fluctuations on model dynamics may be more significant than that of the 
delay. G. Meenakshi and P.R. Maheswari (2020)have derived the posterior distribution using the prior information in the 
Bayesian methodology. In the viral history, replication is increasing nature according to the time. So the viral replication is 
considered as Rayleigh random variable. The predictive distribution obtained and finding the expected number of viral 
replication is illustrated through the graph. Steven G. Deekset al (2021) have developed the success of antiretroviral 
therapy (ART) for people living with HIV, lifelong treatment is required and there is no cure. HIV can integrate in the host 
genome and persist for the life span of the infected cell. G. Meenakshi and P.R. Maheswari (2021)have derived the 
posterior distribution of Weibull distribution by using various priors and estimate the Bayes risk. Simona Claudia Cambrea 
(2022) have modelled Comparing different variables of HIV-positive pregnant women from the two HIV-AIDS CRs, there 
are significant differences between the mean value of hemoglobin, CD} level, environmental area, marital and amniotic 


membranes status, and HIV patient stage in the last trimester of pregnancy (p < 0.05), but without any differences in 
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mother’s mean age, education level, type of delivery, breastfeeding, the duration of cART administration, HIV viral load, 
and survival rate. In this research the seven stage of viral life cycle is considered as the markov differential equation and 


viral replication is estimated by Bayesian methodology. 
3. STOCHASTIC MODEL FOR SEVEN STATES IN THE VIRAL LIFE CYCLE 
1. Binding 


During the first Stage of HIV life cycle, the virus binds to receptors on the surface of CD}T cells. CD{TCells, also call 


helper T cells, are a type of White blood cell that alerts immune cells that there is an infection in our body. 
2. Fusion 


HIV is an enveloped virus, meaning that its genetic information is protected by both a protein shell and a lipid layer called 
an envelope.Once HIV binds to receptors on CD}T cells, it initiates the fusion of its envelope with the membrane of the 
CD{Tcell using a glycoprotein called GP120 Trusted Source. Glycoproteins are molecules made of chains of 


carbohydrates and proteins. Fusing with the membrane of your CD}T cells allows the virus to enter the cell. 
3. Reverse transcription 


Reverse transcription is a process of converting genetic information in the form of RNA into DNA. RNA and DNA contain 
similar genetic information but are structurally different. RNA is typically made up of one long chain of genetic 
information, while DNA is made up of a double strand. The virus converts its RNA into DNA by releasing an enzyme 


called reverse transcriptase. This process allows the virus’ genetic information to enter the nucleus ofCD}T cell 
4. Integration 


Once HIV has converted its RNA into DNA, it then releases another enzyme called integrase inside the nucleus of CD{T 
cell. The virus uses this enzyme to combine its DNA into the DNA ofCD}T cell.At this point, the infection is still 


considered latent and is difficult to detect even with sensitive laboratory tests. 
5. Replication 


Because HIV is now integrated intoCD}T cell’s DNA, it can use that cell’s machinery to generate viral proteins. During 


this time, it can also produce more of its genetic material (RNA). These two things allow it to create more viral particles. 
6. Assembly 


In the assembly stage, new HIV proteins and RNA are sent to the edge of CD} T cell and become immature HIV. These 


viruses are non-infectious in their current form. 
7. Budding 


During the budding stage, the immature viruses push out ofCD}{T cell and then release an enzyme called protease that 


modifies proteins in the virus and creates a mature and infectious virus. 
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In the viral dynamic study, the maturity of the viral replication during the seven stages, namely (1) Binding (2) 
fusion (3) reverse transcription (4) integration (5) transcription (6) translation (7) viral assembly and maturation. Those 
seven stages are considered as continuous time, discrete state Markov chain as follows, 
t;;denote the time duration between the i" stages to j" stage i= 1,2,...;7,] = 2,..,7. 

fiz > ty3 > bq > bas > bs6 > be7 


Let P; j(ti i) denote the probability in the i*” stage (state) to the (j*"state) at the time durationt; j- The probability is 


decreasing nature for every succeeding stateP;; > 0. 
Where 


The Expected time required per new virus maturation is denoted by E(T) = Yi; tijPi;, P,; > 0,andt;; > 


0, i #j The seven stage (state) transition probabilities are as follows 
P, =Pi2 
P, = Py2 Pag 
P; = Py * Pog + Pag 
Py = Piz * P23 * Paq* Pas 
Ps = Py2 * Pez * P3q* Pas * Pse 
Pe = Pz * Pez * P34 * Pas * Ps * Pez 
Where 
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Figure 1. 


The seven state transition probability matrix for nCD}T cell is given by 
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In the viral lifecycle, the Expected time of new virus budding becomes the modified Kolmogorov differential 


equation it is denoted by 


Pio(tiy) = Pa Pi (ti) X Ps,0(te) 
f=} 


1 


i 
s 
is the probability that i? state will be in the O'” state at time ty. + tz3 + ts, + tas + tse tte = t, where tj; are the 


continuous time interval will be decreasing nature of succeeding periods it is given by 
Pig(t) = P(X =) = 0/X(to1) =i, X(t2) = J, XCy3) = k, X(tg4) = L X(t) = m,X(t67) = n) 


This model follows the modified linear growth process with immigration growth rate (when immigration growth 
rate is the dependent to one stage to another stage. Where ¢j;the intervals are not equal. But deceasing random variable 


with interval 
[eg ty] t AO, 1)2 fF, jE cce,7. 
[ao,a,] # [a;, a2] # [az, a3] # [a3, a4] # [ay,a5] # [as,a6] # [a6, a7] 
[ao,a,] > [a;,a,] > [a,,a3] > [a3,a4] > [ay,a5] > [as,a,] > [ag,@,] then its density function is denoted by, 


2et 
a(ite—*)’ 


f(t) = 


a; > 0, a, > Qa, Py mae > az. 


a; —is theevery stage varying ratio (RNA concentration) Expected time to maturation of viralreplication. It is 
given by 
a7 
E(T)= | t-f(@)-dt 
a 


az 


t 2e! 


E(T) = aie 


C 
= 2 (27 **— dtLetz =i+et 


aa, (it+e*) 


2 ¢ logz 
== | A dzlogz =t 


ay 


2 E 
a 


2 
| dz=e'dt 


E(T)= “ (logz)? 


= - [log(1 + e*)] 
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E(T?) = | 2t* 6° , 
~ F a(1 +e) 
1 
_ 2 a7 tel _ t 
= Fle Gren etz =l+e 
az 
2 ¢ 2logz 
=={ = dz2logz = t? 
az 
4 logz 
=-| = dzdz = e‘dt 
ay 
_ 4[(ogz)? 
a 2 
2 
= _ thogz)" Ta; 
== (logz)? 
=o 0gz 


= hlog(t +e')]? 


V(T) = E(T*) — (E(T))? 


2 logz)?) 
= —(logz)* =|") 
a a 
= ieee —— 
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4 4 
V(T) = —log(1 +e") — log(1 +e") 


4. BAYESIAN MODEL FOR VIRAL REPLICATION 


Viral replication in the CD{T cells of Biological system is similar but it is the reverse order to the chemical reactor 
Experiment. Developing viral replication model using the new modified Exponential power distribution. In the chemical 
reactor experiment outlet concentration is smaller than the input concentration. But viral replication system, outlet 
concentration is greater than the inlet viral concentration. In the new virus maturation, there are seven stages. The average 


viral replicationE (y) = n(6;/€) is the non-linear function. 
Where, 
E(y) = n(6;/é) 
y = n(9;/§) + e; 
n(6;/5) = iF 
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y- is the outlet viral replication per CD}T cells after infection. 
€- is the inlet viral fixed viral RNA to the CD}T cells. 
@- is average replication for 7" stage, 
6 = f (G,, Og, .. 07, X4, Xz, X7). 

= O,X, + OX. ++ +O5x7, +€i=1,2,...,7 
Where, 
x; — i stageCD}{T cell DNA and enzymes concentration 
6; — i’ stage viral RNA Concentration, 

Xp = 2x 4,0 = 1,2, ....7 and xq = 2x5; Xgl 
i‘’StageCD}T cell concentration is depends on the provirus stage concentration 
6,andx; are the random variable with parameter 0,0 , 6 
y -Outlet viral replication. 

o- is the scale parameter in each stage variation of enzymes concentration. 
B- is the nuisance parameter, 0 < 8 < 1 possitively correlation between CD} T cells 
DNA Concentration and viral RNA concentration for all the seven stages. 


At the time of infection, viral RNA concentration in the CD{Tas unity. At the next stage concentration may be 
similar to the chemical reactor experiment but it results are in the reverse order. In the chemical reactor experiment, outlet 
concentration is random variable with its range is—oo< y < «, At the time of outlet, the concentration is reduced in the 
chemical reactor Experiment the ever follows the exponential power distribution is used by Diananda (1949), Box (1953b) 


and Turner (1960). They assumed that input concentration is fixed and outlet concentration is random variable. 


In this research, the final stage viral concentration is increasing nature. It is follow the modified exponential 


power distribution with range 0 < y < with parameter 6,¢, B. 


y—E(Q) =e 
The likelihood function is given by 
,,2(B+1) 
1(0,0,B/y) & 1/7 w(B)e"exp {-cw) A 
Where, 
e = yy — EO), O0<Bsl1, O<y<u, o>0. 
[r3/,(4+p)]) "8 
C(B) = 1} I = 
[T/20+B)] 
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y 

r3/,(1 +) °? 
w(B) = 1/, he 2) 3/ 
(1+pyr4/,(1+p)} 


The non-informative the prior distribution of o is given by P(a) = Y/q. 
The posterior distribution of @ given o, B is propositional to 


P(6/o,B,y)a l(y/8,0, B)p(a) 


,2(B +1) 
ao” exp) |= 2 
oO oO 


a nti) 


(6) <> rej\20+B) 

e e 

e > (-) 0<P<1, o>0O, 
t=1 


The Prior distribution of 8 is uniform distribution 
P(B/y) =1 O0<f<l 
The marginal density of the 6 is integrated out (1/ a) 


The posterior density is given by. 


o 


n 
» 4 e248) 
Peo/Byya [ ree S (4) d(1o). 
. i=1 
Let (1/6) =X 
n 
a | xirte-B) Y (ex 0+) dx. 

_ i=1 


The marginal distribution becomes the weibull distribution. 


t 


P(O/B,y)a | xMt2-1 e-C(BYNEX) dy, 
0 


Where a = c(f) (> e;)°,b =n +2,c = 2(1 +8). 


The normalizing constant 


x 


x (t2)-1 e~c(Be(nx) dy. 
0 


~ C(B)OS e)*. (n + 2) 


1 
~ c(B)(E yj —n(0/€)).2(1 + B)(n + 2) 


Impact Factor (JCC): 6.2284 NAAS Rating 3.45 


Stochastic Modelling of Seven Stages Viral Life Cycle and Estimation of Viral Replication Using 59 
Modified Exponential Power Distribution-A Bayesian Approach 


lial te: ey yi-n(8/)(/o) 


P(6/B,y) ~ n(n t+ De(BA(y% —n@/D) 


Mean of the distribution is given by 


Nebaia 


E(X) = C(B)() (1-6/8) (1+ Yn 42) 


6. Numerical Result 


Table 1: Maturation of Viral Replication for Different Time Periods 
Mean and Variance 


_ Time | Mean 
I 1.1406 1.5715 
2 1.8474 5.5416 
3 2.6479 17.4375 
4 3.4901 51.1080 
5 4.3487 144.064 
6 5.2136 398.21 
7 6.0809 1090.5 


E(x) &Vx 


Graph 1: Maturation of Viral Replication Various Scale 
Parameters Mean and Variance 


Graph — I that Maturation of viral replication is mean and variance increases. Whenever the time increase. 


Table 2: Viral Replication for different Time Periods. “",” Vi"'i"us s~’"*"Parameter of the 
Exponential Power Distribution, #: = 7,2, ...,7""* 


0.8 : 1.6 
1 1.2 1.4 1.6 1.8 
1.2 1.4 1.6 1.8 2 
1.4 1.6 1.8 2 22 
1.6 1.8 2 2.2 2.4 
1.8 2 2.2 2.4 2.6 
: : 2 2.2 2.4 2.6 2.8 
Ky) 57.12 75.66 92 137.91 160.05 205.15 232.12 
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Inlet Viral __ 0.32 0.37 0.4 0.54 0.57 0.67 0.7 
concentration(#) 
Fe cceeaisiaaad ia 75.66 | 92 | 137.91 | 160.05 | 205.15 | 232.12 | 
E(y) 
250 
& 
Y 200 
Z 
6 150 4 
% 
3 
5 100 + 
a 
=F 
@ so + 
2 
0 
0.32 0.37 04 0.54 0.57 0.67 0.7 
Inlet viral concentration 


Graph 2: Viral Loads Various Scale Parameter of Exponential 
Power Distribution. 


The viral RNA Concentration increases, the Average viral replication of outlet also increases. 


Table 3: Viral Replication for Different Time P“,""ids""I ‘ie V"“"""’s Shape" .zrameter of the Exponential 
Power Distribution,o« = 1, 2,...,7 csauare Ofxa) 


1 0.3 0.5 0.8 1 1.2 1 
3 0.5 0.8 1 1.2 1.4 1.6 1.8 
9 0.8 1 1.2 1.4 1.6 1.8 2 
81 1 1.2 1.4 1.6 1.8 2 2.2 
6561 1.2 1.4 1.6 1.8 2 2.2 2.4 
43,046, 1.4 1.6 1.8 2 2.2 2.4 2.6 
721 1.6 1.8 2 2.2 2.4 2.6 2.8 
1.85302 
O19E+1 
5 
0.32 2.96483236E+15 
0.37 3.33543641E+15 
0.4 4.07664450E+15 
0.54 5.87357250E+15 
0.57 6.44612046E+15 
0.67 7.56264105E+15 
0.7 8.29543964E+15 
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Average outlet virus E(y) 


4.00E+15 4 
3.50E+15 — 
3.00E+15 — 
2.50E+15 — 
2.00E+15 — 
1.50E+15 + 
1.00E+15 — 
S.00E+14 ~ 


0.00E+00 


0.32 0.37 o4 0.54 0.57 0.67 0.7 


Inlet viral concentration ¢ 


Graph 3: Viral Loads Various Scale Parameter of Exponential 


Power Distribution. 


The inlet viral concentration in terms of square of x; increase.The outlet viral concentration also increases, of the some 


period as compare to the previous table. 


Table 4: Viral Replication for Different Time Periods. , :V r ous{ _Parameter of the 


a = i fa Tse a 
Exponential Power Distribution, e:¢ = 1,2,...,7- 
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Average outlet virus E (y) 
a 
o 


——“ Se 1. 

1 

3 0.5 0.8 1 1.2 1.4 1.6 1.8 

6 0.8 1 12 1.4 1.6 1.8 2 

9 1 1.2 1.4 1.6 1.8 2 Ds 

12 1.2 1.4 1.6 1.8 2 oe) 24 

15 1.4 1.6 1.8 2 Be 24 2.6 

18 1.6 1.8 2 2 2.4 2.6 2.8 

E(y) 25.53 35.26 | 3632 | 64.04 74.89 96.61 109.61 
2 | 032 | 037 | 04 | 054 | 057 | 067 | 07 | 
E(y) | 25.53 | 35.26 | 3632 | 64.04 | 74.89 | 96.61 | 109.61 | 

120 5 


0.32 0.37 o4 0.54 0.57 0.67 0.7 


Inlet viral concentration § 


Graph 4: Viral Loads Various Scale Parameter of Exponential 


Power Distribution. 
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The inlet viral concentration increase, the outlet viral concentration (replication) also increasesafter some period. 


Table 5: Posterior Density of Viral Load Fr%m ¢™C*r''{ R@lication For 
Different Time Periods with Special Caseo: = a er The Scale 
Parameter of the Exponential Power Distrinution 


10 0.129 
20 0.155 
30 0.220 
40 0.329 
50 0.44 

60 0.503 
70 0.791 


Graph 5: Posterior Density of Viral Loads Using Various Scale 
Parameter of Exponential Power Distribution. 


The gradually increases the posterior density of viral replication. 


Table 6: Posterior Density Viral Load for° he’ i 1 Req ication for 
Different Time Periods with Special case « = vi Pa = >l'the Scale 
Parameter of the Exponential Power biseridution 


[ee PT 
10 0.150 
20 0.158 
30 0.178 
40 0.486 
50 0.507 
60 0.524 
70 0.706 
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Graph 6: Posterior Density of Viral Loads Using Various Scale 
Parameter of Exponential Power Distribution. 


The gradually increases the posterior density of viral replication. 
Table 7: Posterior Dei‘ity*”" 4°" for the Viral Replication for Different Time Periods 


with Special Cases = \,iral Load |The Scale Parameter of the Exponential Power 
B@ O75,0=2..... ‘ 
Distribution 


a a 2 
10 0.129 
20 0.169 
30 0.237 
40 0.395 
50 0.480 
60 0.490 
70 0.701 

08 5 
0.7 4 
06 + 


P(6/B.y) 


0.3 4 
0.2 4 
0.1 4 
0 
10 20 30 40 50 60 70 
ei 


Graph 7: Posterior Density Viral Loads Various Scale 
Parameter of Exponential Power Distribution. 


The gradually increases the posterior density of viral replication. 
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Table 8: Mean of the Distribution of Viral Replication For 
Different Time Periods with Special Case N=7 
a Mean of the Distribution 

1 2.47 
2 4.94 
3 7A 
4 9.88 
5 

6 

7 


12.35 
14.82 
17.29 


Mean of the 
distribution 
e 


ei 


Graph 8: Viral Replication for Various Scale Parameters. 


The error term of the distribution is increase. Average viral distribution positively correlated. 


7. CONCLUSION 


In this research, the seven stage of viral replication are modelled as kolmogrov differential equation and the average viral 


replication is considered as non-linear function. Then the posterior distribution of seven stages using Exponential Power 


Distribution and obtained the average distribution of viral replication for the future period through the posterior 


distribution. It is very much useful for the treatment to the HIV infected person. 
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